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Abstract 

In mobilesat applications, frequency offset on the re- 
ceived signal must be estimated and removed prior to 
further modem processing. A straightforward method 
of estimating the carrier frequency offset is to raise the 
received MPSK signal to the M— th power, and then 
estimate the location of the peak spectral component. 
An analysis of the lower signal to noise threshold of 
this method is carried out for BPSK signals. Predicted 
thresholds are compared to simulation results. It is 
shown how the method can be extended to 7 r/M MPSK 
signals. A real-time implementation of frequency offset 
estimation for the Australian mobile satellite system is 
described. 

1 Introduction 

Frequency uncertainty in mobile satellite applications 
arises from doppler offsets due to vehicle motion or local 
oscillator drift. As the symbol rate r 5 is usually low, the 
frequency offset can be a substantial fraction of r s and 
must be removed prior to receive filtering and carrier 
or timing synchronisation. 

A number of previous approaches have used fre- 
quency offset detectors in a feedback loop to correct 
for frequency shifts e.g. [1]. Relatively little attention 
has been given to feedforward estimators which operate 
over a block of symbols. In burst mode applications 
where rapid acquisition is important, the feedforward 
form of estimation is desirable. 

Estimation of the frequency offset during carrier 
only or preamble portions of the received signal is eas- 
ier than when modulation is present. Most of this paper 
concerns the case when random BPSK or QPSK sym- 
bols are transmitted. Estimation of the frequency offset 
with carrier is examined in [2]. 

A relatively simple form of estimating the fre- 
quency offset, involves stripping the modulation by rais- 
ing the MPSK signal to the M — th power, and then es- 


timating the location of the resulting spectral peak via 
the Fast Fourier Transform (FFT) algorithm. This will 
be called the FFT method. Other forms of nonlinear- 
ity may be used for removing the modulation (e.g. [3]), 
but only the M — th power operation is considered in 
this paper. 

The frequency resolution of the FFT method is lim- 
ited by the the duration of the received signal. How- 
ever it is possible to improve the accuracy of locating 
a single spectral component by interpolating between 
the largest bins, or by evaluating the Discrete Fourier 
Transform iteratively to locate the spectral peak. In 
practice, small errors in the frequency estimate (e.g. 
less than 1% of r*) may not matter since they can be 
removed during carrier phase recovery. 

An interesting alternative to the FFT method has 
been reported in [4], This technique effectively esti- 
mates the autocorrelation function of the M — th power 
signal at selected lags and uses this to efficiently esti- 
mate the spectral peak location. 

Section 2 analyses the FFT method for BPSK, in 
order to predict the probability of failure. In section 
3, simulated results are compared to the predicted re- 
sults for BPSK. In addition simulations of QPSK are 
presented. In section 4 a modified version of the FFT 
method is described which can be used with w/Af ver- 
sions of MPSK. Finally, section 5 presents a summary 
of the real-time implementation of frequency estima- 
tion and signal detection strategies for preamble and 
message signals. 

2 Threshold Analysis for the 
FFT Method 

An important question concerning the FFT method of 
frequency offset estimation is when the technique is 
likely to fail due to a noise peak in the FFT being larger 
than the desired signal peak. In this case, since the of- 
fending noise peak may occur at any bin in the FFT 
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Figure 1: Baseband Signal Model for FFT Method of Frequency Offset Estimation 


output, the resulting estimate is effectively a uniformly 
distributed random variable with no correlation to the 
actual offset. The probability of failure is primarily a 
function of the Ei,/Nq , the number of symbols used, 
and the signal constellation size M . This failure thresh- 
old is analysed for BPSK signals in this section. 

The signal model employed is shown in Figure 1. 
Thick lines indicate complex- valued quantities. Stan- 
dard square-root raised-cosine transmit filtering is as- 
sumed, with transmit filter impulse response /*(/). At 
the receiver, frequency estimation must be carried out 
before normal receive filtering since the frequency offset 
may be a substantial fraction of the symbol rate. In this 
model an ideal filter with twice the normal bandwidth 
is used, i.e. f c = r s Hz, prior to frequency estimation. 
It is assumed that the filtered signal is sampled and 
then decimated by a factor of 2, to two samples per 
symbol period T. This is the minimum possible sam- 
pling rate, and furthermore ensures that the noise sam- 
ples are independent which makes the following anal- 
ysis tractable. The use of 2 samples per T, allows a 
frequency offset f 0 of up to ±r s /4, for example, for 
QPSK, or twice this amount for BPSK. The effect of 
not decimating the signal, and using 4 samples per T, is 
investigated by simulation in section 3 and found to be 
quite similar to the minimum sampling-rate case under 
discussion. 

Received signal samples before the squarer can he 
expressed 

r n = -Sn + U’n ( 1 ) 

with the signal component s n 

s n = a k h(r + » 2 T, - kT ) ( 2 ) 

k 

where a* is the k— th transmitted symbol = ±1, T s is 
the initial sampling period = T/4 and r is an arbitrary 


time offset. The signal samples are taken without tim- 
ing recovery or matched filtering and so contain signifi- 
cant InterSvmbol Interference (ISI). The squared signal 
is a “noisy" complex exponential; this section assumes 
that the random components in s 2 are un cor related and 
small compared to the thermal noise terms. These as- 
sumptions appear to he justified at low signal to noise 
ratios. <r| will denote E'{|sn| 2 } = £&/4. The zero-mean 
complex-valued noise samples w„ have real and imagi- 
nary parts whose variance will be denoted equal to 
Aq/T 

Consider now the power estimates from the Dis- 
crete Fourier Transform of A r received samples raised 
to the M—th power i.e. r 2 , rj - ■ - r%_i . Let 

,v-i 2 

P m = £>„ + «■„ IV' 2 "””/" (3) 

n — 0 

for m = 0, • • • A r — i Let the bin containing the signal 
peak be at m — m$; this value of P m will be denoted 
P$. The other bins contain independent noise power 
estimates. Any noise bin (m ^ ms) will be denoted 

/V 

To determine the mean value of the noise power 
estimates, from the last equation 

E{P n ) = £|l^l>P + «’r) 2 K + “p 2 

V P 7 

p _y 2l r m ( ; .-,)/.v j for m ^ ms 

where * denotes complex conjugate. The only non-zero 
terms in this expression arise from one of the signal 
times noise components and the noise-only components. 
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which, remembering that, the noise samples are 
pendent, simplifies to 

£4/?{|s r ne{K,n = 

r 

and the latter term is 
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which simplifies to, using E { J' 1 } = 3<x A r when 

mo 
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Figure 2: Failure Probability for BPSK Signals 


r 


The final result is 

E{Ps} = SN(al(rJ, + (r%) (5) 

Space does not permit a full derivation of the mean 
and variances of P$ and P\ in this paper. Proceeding 
in a similar fashion to that illustrated above, the fol- 
lowing results are obtained: 


°7V 

= E{P*} 

(6) 

E{Ps} 

= A'M + f -{Pn) 

(T) 

<TPs 

= <7p v + l6N 3 ir- H ff 4 s {al + a- R ) 

(8) 


These results assume that terms of order N have 
been neglected from the variance expressions. 

In order to compute the probability that Pjv > Ps 
we need to know the density functions of Ps and P5. 
These distributions are not. obvious due to the non- 
linear operation before the FI T. However it may be 
observed that (5) to (8) correspond exactly to the ex- 
pressions for mean and variance of chi-square and non- 
central chi-square distributions (e.g. [5] (1.1.109) and 
(1.1T 22)) each with 2 degrees of freedom. Using these 
distributions 

PpJpn) = Ee-rN/^ 2 for PN > 0 (9) 

PPs(Ps) = 

for ps > 0 (10) 

where 2 cr ? = E{Ps}, s 2 = A r2 cr$ and /q() is the zero-th 
order modified Bessel function of the first kind. 

Let p F denote the fail probability when any Ps > 
Ps ■ Then 

p F = 1 — prob( all Pv < Ps) 


For a given value ps, of P5, since the Ps are indepen- 
dent, 

prob( all Ps < ps) = (^J Ps[PN)dpN > j 

and so 

p F = 1 - J - e -rs / J ' 7 ) P/v (p.s)dp.v (11) 

When (10) is substituted i Jit o (11), the failure proba- 
bility may be integrated numerically. 

From here it is easy to predict the variance of the 
resulting errors in the frequency estimate f 0 . Letting 
Af 0 — f 0 ~ f o : as A/ 0 is uniformly distributed, with 
either \Af 0 \ < /.</ 2 when the FFT fails, else |A/ 0 | < 
f s /2N, gives 

= ti( ( t )V+ 4 )!<| - w) ) (12) 

where TV is the FFT size. 

3 Simulation Results for the 
FFT Method 

The signal model in Figure 1 was simulated in order to 
check the predicted FFT thresholds derived in sect ion 2. 
A normalised symbol rate of 1 symbol per second was 
used. The continuous-time signals in Figure 1 were 
simulated with 4 samples per T. The excess bandwidth 
parameter in the transmit filter was 40%. 

During the simulations, packets of random bits 
were generated and a random frequency offset between 
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Figure 3: Frequency Accuracy for 61 BPSK Symbols 


Figure I : Frequency Accuracy for 64 QPSK Symbols 


±r,/4 was added. The FFT size was twice the number 
of symbols received, without any zero padding of the 
data samples and with the decimation shown in Figure 
I . When the data was zero padded with an equal length 
of zeros, the FFT was four times the number of symbols. 
The zero padding gives improved performance since it 
reduces the signal loss caused by frequency offsets at 
non-integer multiples of the initial frequency resohit ion, 
and the frequency accuracy is improved. It could be 
argued that the effective number of independent noise 
bins in the FFT output is the same (A r ) with or without 
zero padding; this seems to be supported by simulation. 

Figure 2 shows a plot" of FFT failure probability 
for 32, 64 and 128 bit BPSK signals. The predicted 
values from (11) are plotted with solid lines. These 
simulations used twice oversampling (i.e. 2 samples per 
T) and zero padding. (Since the analysis effectively 
assumes bin-centred doppler frequencies, the use of zero 
padding is appropriate to remove most of the loss due 
to random frequency offsets.) The figure shows fairly 
good agreement between the predicted and measured 
results for BPSK. The lack of simulation results at low 
values of pr reflects the difficulty in simulating rare 
FFT failures at high signal to noise ratios. Plots for 
7t/2 BPSK are discussed in the next section. 

In Figure 3 the standard deviation of the frequency 
error A f 0 is plotted for 64 bit BPSK signals, and com- 
pared to the predicted value from section 2. In the 
evaluation of (12), N was doubled for comparison with 
simulations using zero padding. All the simulations 
show approximately similar thresholds in the vacinity 
of Ei>/Nq = 4 dB, below which the frequency error 
increases rapidly as the failure probability becomes sig- 
nificant. It can be seen that twice oversampling without 
padding gives slightly worse results, as expected. 

Four times oversampling, with padding, has been 


included for comparison. In this case the noise samples 
are correlated and the noise spectrum is no longer flat, 
as it is for section 2. Notice that the final accuracy, 
above threshold, is the same as two samples per T with 
padding, as expected. 

Figure 4 shows 64 symbol QPSK simulations of fre- 
quency accuracy using the FFT method. Again the 
vertical scale is relative to r, — 1 symbol/sec. Tie 
threshold is significantly higher than the M — 2 case 
(almost 6 dB). Again the thresholds for 2/\, or 4r. s sam- 
ple rates appear to bo similar, and the advantages of 
padding can be seen. 

Reliable frequency estimation below the thresholds 
shown in figures 3 and I could be obtained by using 
more symbols per FFT, or by combining results from 
separate FFTs (see section 5). 


4 Frequency Estimation for 
t t/M MPSK 

Variants of MPSK modulation such as 7t/2 BPSK (or 
Aviation BSPK) and t/ 4 QPSK can readily be han- 
dled by the FFT method of frequency offset estimation. 
These modulation schemes are now used in some mo- 
bile satellite systems (e.g. [6]). They rotate alternate 
symbols by n/M radians which reduces the amplitude 
variations in the transmitted signal. 

The effect of these symbol rotations can be appre- 
ciated as follows. Suppose that the receive signal has 
been ideally filtered and sampled at one sample per T 
with no ISI or thermal noise present. The sampled sig- 
nal could be written 




a^eV 4 , o 2 e 24 . ase^V 34 ,--- 
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Figure 5: Sample Power Spectrum for t/2 BPSK 


where <5 = 2 xf 0 T represents the phase change per sym- 
bol period due t.o frequency shift. After the A/ — th 
power operation the samples become 


a o > 


§ M 


J* e jM6 a M e mr* 


If (Ik G e j2?r / A/ , k = 0, 1, • • • M — 1, then = 1 so the 
previous sequence can be written 

\ _ e j A f & c i 2 Af <5^ _ j r j3Af ^ ... 


This represents a discrete-time complex exponential, 
whose frequency is M times the original / 0 , which 
lias been modulated by a sequence of alternating sign, 
---1, — 1, 1, — 1 - of period r 9 / 2. (The MPSK case 
would be the same except for the alternating signs.) 
The signal spectral component is therefore split into 
two discrete components, separated from M f 0 Hz by 
±r\< /2. 

In the pract ical case with more than one sample per 
symbol period and LSI and/or noise, the effect, in the 
frequency domain is similar. Figure 5 shows a typical 
spectral estimate from the FFT method for 128 bits 
of tt/2 BPSK at Eb/No = 6dB. The frequency offset 
was 10% of ?* s ; two components can therefore be seen at 
0.2r,±r*/2. It can be observed that the noise spectrum 
is flat as 2 samples per T were used in this simulation. 

In order to estimate the frequency offset for x/M 
MPSK signals, the peak sum of FFT bins separated by 
r s Hz can be located. This is a straightforward exten- 
sion of the normal peak search method for MPSK. At 
this stage the analysis described in section 2 has not 
been extended to tt/2 BPSK, although this should be 
possible. Some simulation results for tt/A/ MPSK are 
shown in figures 2, 3 and 4. The higher failure proba- 
bility in figure 2 might, be expected due to the spectral 
peak splitting and consequent, lower noise immunity. 


Thi v <r sf,, thresholds are slightly higher for tt/2 BPSK. 
although tins appears 1o hr reversed for x/A QPSK with 
four t imes oversampling. Further simnlat ions would b< 
necessary for precise local ion o( these thresholds since 
a wry small number of FFT failures can significantly 
effect the shape of the frequency error in the threshold 
region. 


5 Real-Time Implementation 

A real time implementation of a frequency estimator 
was required for use in mobilesat(TM) mobile termi- 
nals. The mobile terminals are required to acquire 
phase acquisition of a 3300 symbols/sec tt/ 4 QPSK 
voice activated carrier using a 40ms preamble, or else, 
assuming loss of the preamble, (due to blockages, fad- 
ing etc.), within 180ms. In addition, initial acquisition 
must be obtained using a cont inuous t ransmission (dur- 
ing “Call Set, Up” mode) within 120ms with a 500Ilz 
frequency offset. These acquisition times are specified 
at Eb/No of CdB. Refer to [6] for details. 

Two methods were developed to meet these speci- 
fications. These incorporated detection, frequency esti- 
mation, timing estimation and then phase estimation, 
in that, order. In the case of the preamble, the signal 
was designed to have two discrete components. This 
allowed a very reliable and efficient. PSD correlation 
method to be used for detection, frequency and tim- 
ing estimation. To summarise, four times oversampled 
data was used, and an FFT length of 256 was chosen, 
giving a bin width of 51.6Hz. Detection was based on 
comparison of a peak to average bin correlation ratio 
with a threshold. A detection threshold was chosen by 
experiment. This gave 100% correct detects at values 
of Eb/No down to 3dB. To obtain better resolution 
than the bin width, an interpolation method was used 
between bins. This gave n worst case RMS frequency 
estimate error was of about 6.7Hz, and typically less 
than 3Hz at 3dB. Running on a 50MHz DSP32C, this 
required a processing time of 855ms per estimate, and 
an update rate of every 8 symbols (2.4ms) was chosen. 
This resulted in measured total carrier acquisition times 
(ie including detection, frequency timing and phase re- 
covery) of less than 27ms (typically about 26ms) at 
Eb/No of GdB. About 10ms of this is due to filtering 
and phase recovery delays. 

In the case of the preambleless (i.e. random data) 
signal acquisition, the M = 4 FFT method was used. 
As with the preamide case, four samples per symbol pe- 
riod and a 256 length FFT was used without padding 
(i.e. 64 symbols). This gave an equivalent bin width 
of 12.8Hz (0.0039 relative to ?-,), although using inter- 
polation the frequency resolution could be made much 
less. Frequency errors of 10.3IIz maximum and 3.8IIz 
typically at Eb/No of 3dB were recorded, which were 
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considered accept able given tin- lock range of tin* phase 
recovery technique. 

1 he variance c*f the signal to noise ratio estimate 
was found to be loo high to produce a reliable de- 
tect -ion /reject inn rule when compared with a detec- 
tion threshold. To improve this, the FFT output was 
lowpass filtered to suppress spurious components from 
causing erroneous detections. In addition, hysteresis 
was built into the detection process, to prevent pre- 
mature detection of signals i.e. before the FFT input 
buffer had sufficient, samples of an unfaded signal to 
produce an accurate frequency estimate. These two en- 
hancements allowed a trade-ofT between reliability and 
detection time to be made. This was done experimen- 
tally to meet the mobilesat specifications, and produce 
very reliable detection and estimation. Total process- 
ing delay per frequency estimate was 1050ms and an 
update rate of 8 symbols was chosen. The typical times 
for carrier acquisition were measured to be about 80ms 
at GdB. 

At lower values of E\>fNv . the detection time 
was increased although the 120ms specification was still 
easily met at 3dB. Detection failure was set at about 
1.5dB. Signals with lower Ei,/Nq are too noisy to be 
useful in the given application. 

A crucial performance parameter in the random 
data frequency estimation technique was the input 
noise bandwidth. Due to the high tolerable frequency 
offsets (< 1500Hz) the received signal was digitally fil- 
tered using a low pass filter with cut-off at 1.25 times 
the symbol rate (i.e. slightly larger than in sections 2 to 
4). This was found to give good performance under all 
conditions, although improvements could be attained 
with tighter filtering. 

6 Conclusions 

A method of analysing the lower Ei/Nq threshold for 
the M — th power method of frequency estimation for 
BPSK signals has been presented. Simulations show 
reasonable agreement with the predicted probability of 
the wrong FFT peak being chosen. The FFT method 
can be readily be adapted to tt/M MPSK signals, with 
a small loss in performance for tt/ 2 BPSK. Sample sim- 
ulations for BPSK and QPSK signals give some idea 
of the effects of oversampling and zero padding on 
the FFT threshold. An implementation of the FFT 
method, and related processing, has been tested in real- 
time operation, and optimised on mobile satellite chan- 
nels. This includes using the results of multiple trans- 
forms, rather than one large one, to produce reliable es- 
timation and detection at the required operating point. 
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